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We report the first measurement of the differential cross section for the process 77 — >■ 7777 in 
the kinematic range above the r/rj threshold, 1.096 GeV < W < 3.8 GeV over nearly the entire 
solid angle range, | cos 6^*1 < 0.9 or < 1.0 depending on W, where W and 9* are the energy and 77 
scattering angle, respectively, in the 77 center-of-mass system. The results are based on a 393 fb~^ 
data sample collected with the Belle detector at the KEKB e+e~ collider. In the W range 1.1- 
2.0 GeV/c^ we perform an analysis of resonance amplitudes for various partial waves, and at higher 
energy we compare the energy and the angular dependences of the cross section with predictions of 
theoretical models and extract contributions of the XaJ charmonia. 

PACS numbers: 13.60.Le, 13.66.Bc, 14.40.Be, 14.40.Pq 



I. INTRODUCTION 

Measurements of exclusive hadronic fi- 
nal states in two-photon collisions provide 
valuable information concerning the physics 
of light and heavy-quark resonances, per- 
turbative and non-perturbative QCD and 
hadron-production mechanisms. So far, 
we, the Belle Collaboration, have measured 
the production cross sections for charged- 
pion pairs [lH3|, charged and neutral-kaon 
pairs [3|-[5|, and proton-antiproton pairs [6|. 
We have also analyzed D-meson-pair pro- 
duction and observe a new charmonium 
state identified as the Xc2i2P) 01 . Recently, 
we have examined loJ /ip and fpJ 1^ produc- 
tion and also found charmonium-like struc- 
tures in these final states [8|, Q. 

In addition, we have measured the pro- 
duction cross section for the tt^tt'' and 
7777° final states |lCl| - [l2| . The statistics of 
these measurements are two to three orders 



of magnitude higher than in pre-B-factory 
measurements |l3| . opening a new era in 
studies of two-photon physics. 

In the present study, we report mea- 
surements of the differential cross sections, 
dcr/d\cos9*\, for the process 77 —> 7777 in a 
wide two-photon center-of-mass (cm.) en- 
ergy {W) range from the 7777 mass thresh- 
old 1.096 GeV to 3.8 GeV, and in the cm. 
angular range, | cos 6**1 < 1 (0.9) for W < 
2.0 GeV {W > 2.0 GeV). In this analysis, 
we use the rj ^ "/"/ decay mode only because 
the 77 — >■ 7r"'"7r~7r'^ decay mode has a much 
smaller product of efficiency and branching 
fraction. 

The l'^ J^^ quantum numbers of a meson 
produced by two photons and decaying into 
7777 are restricted to be 0'*"(even)+"'", that is, 
those of /j=even Or Xcj=even mesons. A long- 
standing puzzle in QCD is the existence and 
structure of low mass scalar mesons. In the 
/ = sector, we recently observed a peak- 



ing structure at the /o(980) mass in both 
the 77 — ^ TT'^TT^ and 77 — > tt^tt'^ chan- 
nels P, [l^l ■ Our analysis also suggests the 
existence of another /o meson in the 1.2- 
1.5 GeV region that couples to two pho- 
tons |10| . The significant ss component in 
the 77 meson implies a connection of this re- 
action to the K+K~ H and K°K° d pro- 
cesses. 

At higher energies {W > 2.4 GeV), we 
can invoke a quark model. In leading-order 
calculations, the ratio of the r^Tr" or r/rj 
cross section to that of 7r°7r° is predicted. 
Analyses of energy and angular distribu- 
tions of these cross sections are essential to 
determine properties of the observed res- 
onances and to test the validity of QCD 
based models |l4l4l6l| involving qq produc- 
tion and SU(3) flavor symmetry. It is also 
interesting to compare the behavior of 7777 
production with that of K'^K^ and KgKg, 
which have been measured by the Belle ex- 
periment [sl, Q. The cross section for the 
77 ~^ VV process has not been measured so 
far. 

The organization of this paper is as fol- 
lows. In Sec. II, the experimental appara- 
tus and event selection are described. Sig- 
nal yields and backgrounds are discussed 
in Sec. III. Differential cross sections are 
then extracted in Sec. IV. In Sec. V the 
72(1 270), 72(1525) and other possible reso- 
nances are studied by parameterizing partial 
wave amplitudes. The behavior of differen- 
tial cross sections and W dependence of the 
integrated cross sections at higher energy re- 
gion {W > 2.4 GeV) are compared to QCD 
predictions in Sec. VI. Finally in Sec. VII, a 
summary and conclusion are given. 



II. EXPERIMENTAL APPARATUS 
AND EVENT SELECTION 



Events with all neutral final states are ex- 
tracted from the data collected by the Belle 
experiment. In this section, the Belle de- 
tector and event selection procedure are de- 
scribed. 



A. Experimental apparatus 

A comprehensive description of the Belle 
detector is given elsewhere [17| . We mention 
here only those detector components that 
are essential for the present measurement. 
Charged tracks are reconstructed from hit 
information in the silicon vertex detector 
and the central drift chamber (CDC) lo- 
cated in a uniform 1.5 T solenoidal mag- 
netic field. The detector solenoid is oriented 
along the z axis, which points in the direc- 
tion opposite to that of the positron beam. 
Photon detection and energy measurements 
are performed with a CsI(Tl) electromag- 
netic calorimeter (ECL). 

For this all-neutral final state, we require 
that there be no reconstructed tracks com- 
ing from the vicinity of the nominal collision 
point. Therefore, the CDC is used to veto 
events with charged track(s). The photons 
from a decay of the ry meson are detected 
and their momentum vectors are measured 
by the ECL. The ECL is also used to trig- 
ger signal events. Two kinds of ECL triggers 
are used to select events of interest: the to- 
tal ECL energy deposit in the acceptance 
region used by the trigger (see the next sub- 
section) is greater than 1.15 GeV (the "HiE" 
trigger) , or four or more ECL clusters above 
an energy threshold of 110 MeV in segments 
of the ECL (the "Clst4" trigger). The above 
energy thresholds are determined by study- 
ing the correlations between the two triggers 
in the experimental data. 



B. Experimental data and data 
filtering 

We use a 393 fb~^ data sample accumu- 
lated by the Belle detector at the KEKB 
asymmetric-energy e+e^ collider [iq . For 
an early part of Belle data taking, all neu- 
tral final states were not recorded. Thus, 
this data set is smaller than the hadronic 
data sample available at Belle. 

The data were recorded at several e+e^ 
cm. energy regions summarized in Table HI 
We combine the results from the different 
beam energies, because the e~^e~ cm. en- 
ergy is more than twice our 77 cm. energy 
range for any of the beam energies, and the 
beam-energy dependence of the two-photon 



luminosity function is rather small. We gen- 
erate most of the signal Monte-Carlo (MC) 
events and calculate the two-photon lumi- 
nosity function for 10.58 GeV. We then de- 
rive a correction factor for the other beam 
energies. The correction is less than 0.5% 
over the full range of 77 cm energies consid- 
ered here. The signal MC and the beam en- 
ergy dependences are described in Sec. IV. C. 

The analysis is carried out in the "zero- 
tag" mode, where neither the recoil elec- 
tron nor positron are detected. We restrict 
the virtuality of the incident photons to be 
small by imposing a strict requirement on 
the transverse-momentum balance with re- 
spect to the beam axis for the final-state 
hadronic system. 

The filtering procedure ("Neutral Skim") 
used for this analysis is the same as the 
one used for tt^tt^ and t]tt^ studies |lCl| - tl2| . 
The important requirements in this filter are 
the following: there are no tracks originat- 
ing in the beam collision region and hav- 
ing a transverse momentum greater than 
0.1 GeV/c in the laboratory frame; two 
or more photons that satisfy a specified 
energy or transverse-momentum criterion; 
this requirement is satisfied when there are 
three or more photons each with an en- 
ergy above 100 MeV. The performance of 
the ECL tri gger s is studied in detail using 
7r°7r° events [10|- We also study the trigger 
thresholds using the rjri signal samples. 



C. Event selection 

From the Neutral Skim event sample, we 
select 77 — >■ 7777 candidates that satisfy the 
following conditions: 

(1) the total energy deposit in ECL is less 
than 5.7 GeV; 

(2) each photon candidate is required to 
have an energy of at least 100 MeV, 
and events with four such photons are 
selected; 

(3) the event is triggered by either the 
ECL trigger HiE or Clst4; 

(4) either the sum of the energies of 
the photons in the acceptance region 
used by the trigger is larger than 
1.25 GeV, or all four selected pho- 
tons are within this region, where the 
trigger acceptance is the polar-angle 



range, -0.6255 < cos6' < +0.9563, in 
the laboratory frame; 

(5) of the three possible combinations 
that can be constructed from the four 
photons, there is one in which each in- 
variant mass of the two photon pairs 
satisfies 0.52 GeV/c^ < M^^ i < 
0.57 GeV/c^, where i = 1, 2 is an in- 
dex of the two-photon pairs; 

(6) there is no neutral pion combination 
that is constructed from any two of the 
four photons with a x^ smaller than 9 
in the mass-constrained fit; 

(7) the transverse momentum for the 7777 
system |Sp(| is required to be less 
than 0.05 GeV/c. 

A small fraction of events contain multiple 
combinations of the four photons that sat- 
isfy criterion (5). In those events, we take 
only one combination whose residual for the 
nominal 77 mass (771,, = 0.5478 GeV/c^), 



{M. 



772 ^ 



is the small- 



{M^yl - 777^ 

est. 

We then scale the energy of the two pho- 
tons with a factor that is the ratio of the 
nominal 77 mass to the reconstructed mass, 
niri/My-y i. This is equivalent to an approx- 
imate IC (one constraint) mass constraint 
fit in which the relative energy resolution 
(AE/E) is independent of E and the res- 
olution in the angle measurement is much 
better than that of the energy. This is a 
good approximation for the 77's in this mo- 
mentum range. Using the corrected four- 
momenta of the 77 mesons, we calculate the 
invariant mass {W) and the transverse mo- 
mentum (|Epj|) in the e~^e~ cm. frame 
for the 7777 system and apply cut (7) above. 
We select 31655 candidates in the region 
W < 4.0 GeV. 

We define the cm. scattering angle, 9*, 
as the scattering angle of the 77. The e+e~ 
direction is used to approximate the axis for 
the polar angle calculation because the ex- 
act 77 axis is unknown for untagged events. 
The two-dimensional {W, \ cos(?*|) distribu- 
tion of selected events is shown in Fig. [TJ 

The probability for a signal 7777 event to 
have multiple combinations is sizable only 
near the threshold (about 6% at W ^ 
1.11 GeV), but it is small (less than 2%) 
above W > 1.12 GeV, according to the sig- 
nal MC samples. For different choices of 7- 
pair combinations in an event, the W val- 



TABLE I: Data sample: integrated luminosities and energies 



e e cm. energy 


Integrated luminosity Runs 


(GeV) 


(fb-i) 




10.58 


286 


T(45) 


10.52 


33 


continuum 


9.43 - 9.46 


7.3 


near T(IS') 


9.99 - 10.03 


6.7 


near T(2S') 


10.32 - 10.36 


3.2 


near T(3S) 


10.83 - 11.02 


58 


near T(5S') 


Total 


393 





ues are nearly the same, but | cos 9* \ can be 
different. As the angular distribution is ob- 
served to be flat near the threshold, which 
is also theoretically expected, the effect of 
an incorrect choice is negligibly small. 



III. YIELDS OF THE SIGNAL AND 
BACKGROUNDS 

In this section, backgrounds are identified 
and subtracted and the extraction of the sig- 
nal yield is discussed. 



A. Determination of non-?7r7 
background 

There are two kinds of background pro- 
cesses for the 77 — >■ 7777 signal process: non- 
•qri and ^^rjX backgrounds. The non-?/?/ back- 
ground does not contain an ry pair in the 
final state, while the 7777X background in- 
cludes extra particle(s) in the final state in 
addition to the 7777 combination. In this 
measurement, the non-7777 contribution, aris- 
ing from beam backgrounds or other physics 
processes, is the dominant background in 
the final sample. 

We first determine the number of the non- 
7777 background events using the yields of 
7777 mass sidebands. After subtracting this 
background contribution, we check the pt- 
balance distribution for the remaining com- 
ponent; the signal component peaks near 
|Sp*| = 0, while we expect that the rj-qX 
background does not. 



1. Defining the rjr]-mass sidebands 

The 7777-mass sidebands are defined by 
displacing the central points of the mass 
intervals in selection criterion (5) by 
±0.07 GeV/c^. Two kinds of sidebands 
are defined: Sideband A and Sideband 
B. In Sideband A, the central points for 
the two-dimensional mass cut for (M-y^i, 
My-y2) are (here, we assume M-y^i < My^2) 
(0.545, 0.615) and (0.475, 0.545) in units 
of GeV/c^ and the width of the range is 
±0.025 GeV/c2. Sideband B has central 
points (0.475, 0.475), (0.475, 0.615) and 
(0.615, 0.615). When there are two or more 
choices of 7-pair combinations in an event 
that fall in the same sideband box, we take 
the one that is closest to the nominal cen- 
tral point of each sideband box, (771,, — 0.07, 
77T,,;) or (ttt,^, 777.-;, + 0.07) for Sideband A, and 
{rn,j + / - 0.07, 777^ + / - 0.07) for Sideband 
B. This is similar to the multiple candidate 
selection applied for the signal candidates. 
The Al-yy distributions near the signal and 
sideband regions are shown in Fig. [21 

We also calculate W and | cos 9* \ for the 
Sideband A and B candidates by scaling 
Mj^ i to 777,,j (not to rriri + / — 0.07, which 
would change the threshold mass). 



2. Sideband subtraction 

We subtract the sideband yield to obtain 
the signal component with the following for- 
mula: 

y = y, -o.5yM + o.25nB, 




FIG. 1: Two-dimensional W and | cos 5*1 distribution for the rjrj candidates in data. The same 
distribution is viewed from two different directions. The backgrounds are not subtracted. 



where Y is the signal yield after sideband 
subtraction, Yg is the yield in the observed 
events in the signal region, and YbA ( YbB ) 
the yield of the Sideband A (B) region. Here 
we model the non-?y backgrounds with a lin- 
ear distribution in M^-y, for backgrounds 
with both r] non-77 and non-77 non-77 com- 
binations. The possibility of a non-linear 
background component is included in the 
systematic error (see Sec. IV. E). The yield 
in the signal and sideband regions (before 
the sideband subtraction) is shown in Fig.|31 
To obtain the differential cross sections, we 
subtract bin-by-bin in each two-dimensional 
bin of {W, I cos 6** I) with bin widths AW = 
20 MeV and A| cos 61* | = 0.1. Two or five W 
bins are combined later in the determination 
of the final cross sections. Signal leakage 
into the sideband regions, which amounts 
to 2-5% of the signal size Y and is larger at 
small W, is expected according to signal MC 
simulations. This effect is corrected in the 
derivation of the differential cross sections 
by reducing the efficiency. 



B. pt-unbalanced component 

We expect that the background remain- 
ing after the 7/7/ sideband subtraction is very 
small. To confirm this, we examine the 
W dependence of the yield ratio of the pt- 
unbalanced component R defined as: 



R 



y(0.15 GcV/c < I EPtI < 0-20 GeV/c) 



y(IEPt1<0.05GeV/c) 



(1) 



where Y is the yield after the sideband sub- 
traction in the specified pt region. R is plot- 
ted as a function of W in Fig. [H where 
any excess over the signal MC would indi- 
cate a contribution from 7777X background. 
There is such a small excess just above mass 
threshold. We include the effect from this 
possible background source into the correc- 
tion and the systematic error. Non-7777 back- 
ground is much larger near the W thresh- 
old, and this excess may be due to an im- 
perfect sideband subtraction. We apply a 
—3% correction for the background from 
this source for W < 1.2 GeV. In the W re- 
gion between 1.2 GeV and 3.3 GeV, the R 
value is consistent with the signal MC sim- 
ulation. The reason why the experimental 
data seems to be slightly below the MC for 
R in the 1.4 - 2.0 GeV range is not known, 
but the difference translated to the back- 
ground ratio is negligibly small, less than 
1%. For W > 3.3 GeV, there could be much 
larger rjrjX backgrounds. As described in 
Sec. IV. C, we do not report cross section re- 
sults for W > 3.3 GeV, and apply a —3% 
correction for 3.2 GeV< W < 3.3 GeV. 

We conclude that this kind of background 
is less than 2% throughout the W region. 
1.2 - 3.2 GeV; and assign 2% as the sys- 
tematic error for this source for the en- 
tire W region, 1.096 - 3.3 GeV. These fac- 
tors are obtained by assuming a quasi-linear 
I ^ Pt I dependence of the background and 
extracting its leakage into the signal region 
i\J2Pt\ < 0-05 GcV/c), which is approxi- 
mately 1/6 of the yield in the 0.15 GeV/c < 
\J2Pt\< 0-20 GeV/c region. 
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FIG. 2: (a) Two-dimensional Mr^~f\ vs. M^^-g distribution of the four-photon candidates in data. We 
take M^T-i < M~,~fi. The pt-balance cut with pt < 0.1 GeV/c is apphed using the photon momenta 
before the m^ mass correction, which reduces backgrounds. Red, blue and green boxes show the 
signal, Sideband A and Sideband B regions, respectively, (b) A one-dimensional projection of the 
same distribution where the 77 pair on the opposite side is required to be in the signal mass region, 
0.52-0.57 GeV/c^. The vertical red and blue lines show the signal and two Sideband A regions, 
respectively. Note that there are two entries per event in the signal region. 




W (GeV) 

FIG. 3: W distributions of the yields in the sig- 
nal region (closed circles with error bars, Ys), 
and estimates of backgrounds of the rj non-ry 
component (solid histogram, 0.5{YbA — i^6s)) 
and the non-r; non-77 component (dashed his- 
togram, 0.25Yi,_B). Points with error bars 
show the estimated total background, 0.5yi,A — 
0.25ns. 




2 3 

W (GeV) 

FIG. 4: Energy dependence of R defined in 
Eq. ([T} . It indicates the level of rjijX background 
contamination, for the experimental data (af- 
ter the sideband subtraction, points with error 
bars) and signal MC (short horizontal bars). 



A. Effect of e"'"e beam energy 



IV. DERIVING DIFFERENTIAL 
CROSS SECTIONS 



In this section, we present the procedure 
to derive differential cross sections. 



We generate standard MC events for an 
e+e" cm. energy ^/s = 10.58 GeV. 
We compare the products of the luminos- 
ity function and efficiency (L-yy{W)e in 
Eq. (2)) at three different cm. energies, 
9.46 GeV(T(lS')), 10.58 GeV(T(4S')) and 
10.87 GeV(T(5S')) using MC samples. We 
conclude that, taking into account the inte- 
grated luminosities of tlie different cm. en- 



ergies, the correction factors for the lower 
and higher energy samples cancel almost ex- 
actly. Applying the MC results for 10.58 
GeV to all samples leads to negligibly small 
effects of less than 0.5%. 



B. Invariant mass resolution 

We estimate the invariant mass resolution 
of the 7777 system using the signal MC sim- 
ulation. Since we apply an energy rescal- 
ing using the 77 mass, the W resolution is 
better than that for a pure energy measure- 
ment. We find that the invariant mass res- 
olution is about 0.6% near the threshold, 
W = 1.1-1.5 GeV and approaches 1.0% for 
higher W. We confirm that the experimen- 
tal resolution is at most 10% larger than 
the MC resolution from measurements of 
Pt balance in tt^tt^ production and in the 
77' peak in 77 — > 77' — )■ 77. The resolu- 
tion is much smaller than the W bin widths: 
AW = 0.04 GeV or 0.1 GeV. Since statistics 
are low. we do not unfold our results as in 
previous measurements |10| - [l3| . 



C. Determination of tlie efficiency 

The signal MC simulations for e+e^ -^ 
e~^e~rm are generated using the TREPS 
code [13 and are used for the efficiency cal- 
culation at 32 fixed W points between 1.1 
and 4.0 GeV and isotropically in |cos0*|. 
We evaluate the efficiencies separately in 
I cos 0*1 bins with a width of 0.05, and thus 
the angular distribution at the generator 
level does not play a role in the efficiency 
determination. 

The Qmax parameter that gives a maxi- 
mum virtuality of the incident photons is set 
to 1.0 GeV^, while the cross sections for vir- 
tual photon collisions include a form factor, 

fT^^(o,g2) ^<7^^{o,o)/{i + QywY. Om- 

analysis is not sensitive to the form factor 
assumption, since our stringent pt-balance 
cut {\Y.p;\ < 0.05 GeV/c) implies Q^/W^ 
is much smaller than unity; an approximate 
relation Q^ ^ I X]Pt P holds when only one 
incident photon is treated as moderately vir- 
tual and the scattering angle of an electron 
(or a positron) that has emitted the virtual 
photon is small. Using signal MC Simula- 




□ .9 3.5 3.25 



W(GeV) 



FIG. 5: Two-dimensional dependence of the 
efficiency on {W, j cos^*|). 



tion and replacing the Q^/W^ term by ei- 
ther Q^ /rri^ or omitting it entirely, we con- 
firm that the effect of the form factor choice 
on the cross section is less than 0.5%, where 
nip is the p meson mass. 

Samples of 400,000 events are generated 
at each W point and are passed through 
the detector and trigger simulations. The 
obtained efficiencies are fitted to a two- 
dimensional function of {W, |cos0*|) with 
an empirical functional form. 

We embed background hit patterns from 
random trigger data into MC events. We 
find that different samples of background 
hits give small variations in the selection ef- 
ficiency determination. A M^-dependent er- 
ror in the efficiency, 3 - 4%, arises from the 
uncertainty in this effect. Figure [S] shows 
the two-dimensional dependence of the effi- 
ciency on ( W, I cos6'*|) after the smoothing 
fit. 



D. Derivation of differential cross 
sections 

The differential cross section for each {W, 
I cos6'*|) point is given by: 



da 



AY 



d\ cos 6* I AM^ A I cos 6* \ J Cdt L^^ (W) e B^' 

(2) 
where Ay is the signal yield after the 
77-mass sideband subtraction, AW and 
AjcosS*! are the bin widths, J Cdt and 
Lj-y(W) are the integrated luminosity and 
two-photon luminosity function calculated 
with TREPS [lil, respectively, e is the ef- 
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TABLE II: Bin sizes. The lowest bound of the 
first W range (1.0957 GeV) corresponds to the 
mass threshold. 



W range 


AW 


A| cos 61*1 


maximum 


(GeV) 


(GeV) 




Icosri 


1.0957 - 1.12 


0.0243 


0.1 


1.0 


1.12 -2.0 


0.04 


0.1 


1.0 


2.0 - 2.4 


0.04 


0.1 


0.9 


2.4 - 3.3 


0.10 


0.1 


0.9 



ficiency, and B^ is the squared branching 
fraction for 77 — > 77. The over-subtraction 
of signal in the sideband due to the leak- 
age of the signal into the sideband region is 
evaluated in the MC, separately, and finally 
included in the efficiency e. 

The bin sizes AW and A| cos 6**1 and the 
maximum | cos 9* \ for which we obtain the 
differential cross section are summarized in 
Table HT] We first derive the differential 
cross sections for bin widths of AW = 
0.02 GeV and A|cos6'*| = 0.1, and aver- 
age the differential cross section over two or 
five different W regions to obtain results for 
AW = 0.04 GeV or 0.10 GeV, respectively. 

We do not give a cross section for W > 
3.3 GeV. In the W range 3.3-3.6 GeV, 
the charmonium component dominates the 
yield, and we cannot subtract it in a model- 
independent way. We also cannot give the 
cross section including the charmonium con- 
tribution in these bins, because leakages 
from the narrow Xco peak around 3.41 GeV 
into adjacent bins due to energy resolution 
complicate the extraction of cross sections 
in each bin. Above W > 3.6 GeV, we do 
not find any significant signal after consid- 
eration of the backgrounds. 

Figure [5] shows the angular dependence 
of the differential cross sections for selected 
W bins. Figure [7] shows the cross section 
integrated over |cos6'*| < 0.9 for the entire 
W range and that for |cos6'*| < 1.0 in the 
range W < 2.0 GeV. 



E. Systematic errors 

Various sources of systematic uncertain- 
ties assigned for the 7777 signal yield, effi- 



ciency and the cross section evaluation are 
described in detail below and summarized 
in Table Uni 

(1) Trigger efficiency: The systematic er- 
ror due to uncertainty in the thresh- 
old for the Clst4 trigger (- 110 MeV) 
is very small, because photons from rj 
decays have high enough energy. How- 
ever, the efficiency of the HiE trigger 
dominates that of Clst4 except in the 
lowest W region, because the former 
has a looser condition for the number 
of clusters in the acceptance region of 
the trigger. We estimate the uncer- 
tainty in the efficiency for the HiE trig- 
ger to be 4% over the whole W region, 
and treat it as the combined system- 
atic error for the two kinds of triggers. 

(2) ?7 selection efficiency: We assign 6% 
for the selection of the two 77's. This 
corresponds to a 3% uncertainty for 
the efficiency of each 77 reconstruction. 

(3) Overlapping hits from beam back- 
ground and related effects: We assign 
a 4% (3%) error for I^ < 1.44 GeV 
{W > 1.44 GeV) for uncertainties of 
the inefficiency in event selection due 
to beam-background photons, which 
affect the photon multiplicity and rj 
reconstruction. The uncertainty is 
estimated by comparing efficiencies 
among different experimental periods 
and background conditions. We adopt 
the average efficiency from different 
background files, and the uncertainty 
in the average, obtained from the vari- 
ation of experimental yield in different 
run periods, is assigned as the error. 

(4) pt-balance cut: A 3% uncertainty is 
assigned. The pt-balance distribution 
for the signal is well reproduced by 
MC so that the efficiency is correct to 
within this error. 

(5) Sideband background subtraction: 1/3 
of the size of the subtracted compo- 
nent is assigned to this source for each 
bin. We conservatively assign this er- 
ror because we ignore the non-linear 
behavior of the background in the Mj-y 
distribution in the sideband subtrac- 
tion. This effect is expected to be 
large but cannot be determined pre- 
cisely in the lowest W bins. 



10 



J3 

c 

(X> 

w 
o 
u 

:;o 

T3 




\i 






w 


= 1.26 QeV 


i 


H 


*, 












•n 










f 



■\ 


1 


i 


i 


w = 


1.42 GeV 

- I 


1 1 1 



, i i W= 1.54 GeV 

i 





W = 3.05Ge\ 




" 




. ,.H 




• • 1 



|cos e*| 

FIG. 6: Angular dependence of the differential cross sections for nine selected W bins indicated. 
The bin sizes are summarized in Table HIl 
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FIG. 7: (a) The cross section integrated over 
I cos 61* I < 0.9 and (b) over |cos6»*| < 1.0 for 
W < 2.0 GeV. Errors are statistical only. The 
dotted curve shows the size of the systematic 
uncertainty. 



(6) pt-unbalanced background: We have 
applied a —3% correction for this 
background source only in the lowest 
and highest W regions, W < 1.2 GeV 
and W^ > 3.2 GeV. We do not find any 
evidence of such a component, and no 
correction is applied for this effect in 
the other energies. We assign a 2% er- 
ror from this source for the entire W 
region. 

(7) Luminosity function: We assign 4% 
(5%) for W below (above) 3.0 GeV; 
this includes the uncertainties in the 
equivalent photon approximation (3% 
(4%)), the radiative corrections that 
were neglected (1-2%) and the inte- 
grated luminosity (1.4%). 

(8) No unfolding: 

Uncertainty from smearing effects is 
estimated by smearing a modeled res- 
onance function with the W resolu- 
tion and examining apparent changes 
of the cross section. The changes are 
large (^ 7%) only near the slopes of 
the narrowest resonant structure, in 
the region 1.44 GeV <W < 1.60 GeV, 
and smaller (4%) in other W ranges. 

(9) Other efficiency errors: An error of 
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4% is assigned for uncertainties in the 
efficiency determination based on MC 
including the smoothing procedure. 

The total systematic error is obtained by 
adding all the sources in quadrature and is 
11-12% for the intermediate and high W 
regions. It becomes more than 20% for 
W < 1.24 GeV. 

In the resonance analyses for W < 
2.0 GeV in Sec. V, we treat the systematic 
error sources except for (9) as uncertainties 
in the overall normalization, which are cor- 
related in the different (W, \cos9*\) bins. 
For the analysis of the W dependence in 
the high energy region (Sec. VLB), we also 
take into account energy-dependent devia- 
tions for sources (6) and (8). 



V. STUDY OF RESONANCES 

In the total cross section (Fig. [7]), clear 
peaks due to the /2(1270) and /:^(1525) 
are visible along with other possible reso- 
nances. In this section, we first present con- 
sistency checks with previous measurements 
and report improved measurements of some 
of these resonances. 



A. Differential Cross Sections in 
Partial Waves 

In the energy region W < 3 GeV, J > 4 
partial waves (next is J = 6) may be ne- 
glected so that only S, D and G waves are 
considered. The differential cross section 
can be expressed as: 



\Yj- 1 's are not independent of each other par- 
tial waves cannot be separated from the in- 
formation in the differential cross sections 
alone. 

We rewrite Eq. ^ as 



da 



AT:d\ cos 



-(77^7777) 



Y,r 



+Gi\Y,r 



-Gl\Yi\' 



Di\yi? 



(4) 



The amplitudes 5^ Dl, Dl, G§ and Gj 
can be expressed in terms of S, Dq, D2, Go 
and G2 [ll|. Since the square of spherical 
harmonics are independent of each other, 
we can fit differential cross sections to ob- 
tain 5-2, D^, Dl Gl and Gl in each W 
bin. Since \Y2\^ and \Y^\^ are nearly equal 
for I cos 6** I < 0.7 we also fit Gq + G| and 
Gq — G\. Two types of fits are made: the 
"SD" fit and "SDG" fit. G waves are ne- 
glected in the SD fit. 

The spectra of iS^, D^ and D2 obtained 
for the SD fit and G^, G^ and Gg ± G^ for 
the SDG fit are shown in Figs. [5] and HI The 
spectra of S"^ , Dq and Z)| for the SDG fit are 
omitted because they are nearly the same as 
those for the SD fit with somewhat larger 
statistical errors. It appears that the Dq and 
G waves are small enough to be neglected 
in the region of interest {W < 2.0 GeV). 
In that case, S^ and ZJ^ become \S\'^ and 
|D2p! respectively, which simplifies the pa- 
rameterization. In the fits performed here, 
we neglect the G waves completely, and take 
Dq = in the nominal fit. 



^(77 ^W) 



\SYo" + DoY2" + GoY^\^ 
\D2Yi + G2Yi\\ (3) 



where Dq and Go (-D2 and G2) denote the 
helicity (2) components of the D and G 
waves, respectively^, and Yf are the spher- 
ical harmonics in which the helicity A is 
quantized along the 77 axis. Since the 



^ We denote individual partial waves by roman let- 
ters and parameterized waves by italic. 



B. Fitting Partial Wave Amplitudes 

In this subsection, we describe the ex- 
traction of resonant substructure by fitting 
differential cross sections by parameterizing 
partial wave amplitudes in terms of reso- 
nances and smooth "backgrounds". Note 
that we do not fit 5^, D^ and IJ|, but in- 
stead fit the differential cross sections di- 
rectly. Once the functional forms of am- 
plitudes are assumed, we can use Eq. ([3]) 
to fit differential cross sections. We then 
do not have to worry about the correla- 
tions between 5^, Dl and i)|. The S'^, Dl 
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TABLE III; Systematic errors for the differential cross sections. Ranges of errors are shown when 
they depend on W. 



Source 



Error (%) 



Trigger efficiency 

ry-pair reconstruction efficiency 

Overlapping hits from beam background etc. 

Pt-balance cut 

Sideband background subtraction 



4 
6 
3-4 
3 
2-27 (for W>1.2 GeV) 
28 - 60 (for W<1.2 GeV) 
Pt-unbalanced background subtraction 2 

Luminosity function and integrated luminosity 4-5 

Unfolding 4-7 

Other efficiency errors 4 



Overall 



11 - 29 (for W > 1.2 GeV) 
30 - 61 (for W < 1.2 GeV) 



and D2 spectra are compared with the re- 
sults of parameterization. Here we neglect 
the Do and G waves in the fitting region, 
W < 2.00 GeV. 



Quite a few resonances are listed in 
Ref. [23 (PDG) that are known to de- 
cay into 7777 with measured or unknown 
branching fractions to two photons. Be- 
sides the /2(127G) and /2(1525), there 
are /2(1565), /2(1910) and /2(1950) tensor 
mesons, /o(1370), /o(1500), /o(1710) and 
/o(2020) scalar mesons, and spin-4 /4(2050) 
states. So far quantitative measurements 
of the branching fraction to 7777 based 
on observed enhancements in mass spec- 
tra are available for /2(1270) [Ml, [H] and 
/2(1525) J2^. In addition, a phenomenolog- 
ical derivation of the 7777 branching fraction 
based on a K-matrix a ppr oach |24| has been 
tried for the /2(1525) Q. 



To investigate this complicated region, we 
divide our analysis into two parts. First, 
we try to confirm or improve the parame- 
ter T^jB{r]r]) of the well established tensor 
mesons, /2(1270) and /2(1525), by fitting in 
the region W < 1.64 GeV. We then investi- 
gate the higher mass region by fixing most 
of the parameters in the fit from results in 
the low mass region. 



1. Low mass region, 1.12 - 1.64 GeV 

We concentrate on the resonances, 
/2(1270) and /a (1525) by fitting the region 
W < 1.64 GeV. The resonances taken into 
account are the /2(1270), /2(1525) and 
"foiXy, where "/o(^)" is .just a parame- 
terization motivated by the /o(1370) and 
/o(1500). We parameterize partial waves as 
follows: 

S = Af^^Y)e'*- + Bse'^% 
Dq = Bjjo, 

D2 - %(1270)e"^^+^/^(1525)e"^-^+Sl?2, 

(5) 

where A/^(y), Af^(^i270) and ^^^(1525) are the 
amplitudes of the corresponding resonances; 
Bs, B]jQ and B]j2 are "background" ampli- 
tudes for iS*, Dq and D2 waves; (J)y, 02 and (j)^ 
are the phases of resonances relative to back- 
ground amplitudes; (j)^ is the relative phase 
between S and Dq. We set Bjjq = (and 
then (f)s — 0) for simplicity in the nominal 
fit, but we later consider a non-zero Dq con- 
tribution to determine the systematic errors 
for the obtained resonance parameters and 
leave the Bjjq symbol here. 

To parameterize resonances, we use a rel- 
ativistic Breit-Wigner amplitude Aji(W) for 
each spin- J resonance R of mass mji given 

by 
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W (GeV) 



TV" 
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T 



1^5 1^ 1.75 2 2J5 2.5 2.75 3 325 
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+ 



^ttttt 
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FIG. 8: Spectra of S'^ Dq and DJ for the SD 
fit. Those for the SDG fit are nearly identical 
with larger statistical errors. The error bars 
shown are statistical errors that do not include 
correlations. 
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FIG. 9: Spectra of Gg and G^ and Gf) ± Gl 
The error bars shown are statistical errors that 
do not include correlations. 



width Txx{W) is parameterized as [2£ 



AUw) ^ f-^S£±D^ 



ml-W^-imRTtotiW) 



(6) 



For scalar mesons, partial and to- 
tal widths do not depend on VF, while 
for tensor mesons (the /2(1270) and 
/^(1525), /2(1810) and /2(1950)), the 
energy-dependent total width rtot(W^) is 
given by 

rtot{w)^J2^x^iw), (7) 

X 

where X is a tt, K^ rj, 7, etc. The partial 



T^xxiW) = TrB{R^XX) 
D2 {qx{W^)rR) 

V - — 

D2 {qx{m%)rR) 



fqx(W^) 



\qximl) 



(8) 



where r^ is the total width at the resonance 
mass, qx{W^) = y^W^/A^mj^, D2{x) = 
1/(9-1- Sx^ + 2;"'), and tr is an effective in- 
teraction radius that varies from 1 GeV~^ to 
7 GeV~^ in different hadronic reactions [291. 
We assume the same tr value from Ref. Q 
for the /2(1270) and /2(1525). 

For the Air and other decay modes, 
r4^(l^) = Tr,B{R ^ 47r)^ is used in- 
stead of Eq. © for the /2(1270). Param- 
eters of the /2(1270) and /2(1525) are sum- 
marized in Table IIVI T he resonance param- 
eters given in Ref. [20| for the /o(1370) and 
/o(1500) are summarized in Table FVl Back- 
ground amplitudes are parameterized as fol- 
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lows. 



Bs - l3{bs{W-Wo) + cs), 
Bdo = l3HboiW-Wo)+co), 
Bd2 = P''(h2{W-Wo) 



C2 



where /3 is the velocity of the rj meson in 
the c.m.s. and Wo = 2to^. We set Boo = 0, 
that is, &o = Co = 0, in the nominal fit. 
We assume the background amplitudes for S 
and D2 to be real and linear in W to reduce 
the number of parameters. Furthermore, we 
fix arbitrary phases by choosing cs > 0, and 
C2 >0. 

We fit the energy region of 1.12 GeV < 
W < 1.64 GeV. In the fit, we fix the val- 
ues of the parameters of the /2(1270) and 
/^(1525) to those in the PDG [23 except 
for the product T^-yB{riri) for the /2(1270). 

Two hundred sets of randomly generated 
initial parameters are prepared and fits are 
performed for each study. A unique solution 
is obtained with a fit quality of x^/nd/ = 
137.1/119, where ndf is the number of de- 
grees of freedom in the fit. A fit with- 
out f[)(Y) gives a poor fit with ^ /ndf = 
209.7/123. The parameters obtained from 
these two fits are summarized in Table IVIl 
The product r^^6(ry7y) for the /2(1270) is 
ll.S^'^jQ ^ ^i^'l i^ consistent with 12.1 ± 
2.8 eV'in PDG pQl- Figures [TO] to [H] show 
results of the nominal fit to differential cross 
sections, the total cross section, and spectra 
of 52, Dl andL>2. 

Fits where the value of the product 
T^-^Blrj-q) of the /2(1525) is floated while 
that of the /2(1270) is fixed to the PDG 
value, yields three solutions listed in Ta- 
ble [Vn] Thus we fix the former to the PDG 
values in further studies. 

The following sources of systematic errors 
on the parameters are considered: depen- 
dence on the fitted region, normalization 
errors of the differential cross sections, as- 
sumptions on the background amplitudes, 
and the measurement errors of the /2(1270) 
and /^(1525). 

For each study, a fit is made allowing all 
the parameters to float; the differences of 
the fitted parameters from the nominal val- 
ues are quoted as systematic errors. Here 
too, two hundred sets of randomly gener- 
ated initial parameters are prepared for each 
study and fitted to search for the true min- 
imum and for possible multiple solutions. 



Unique solutions are found many times. 
Once a solution is found, several more it- 
erations of the fitting procedure are made 
to confirm the convergence. 

The resulting systematic errors are sum- 
marized in Table IVIIIl Two fitting re- 
gions are tried: one region that is shifted 
lower by one bin (1.08 GeV < W < 
1.60 GeV) and another shifted higher by one 
bin (1.16 GeV <W< 1.68 GeV). Stud- 
ies on normalization are divided into those 
from uncertainties of the overall normaliza- 
tion and those from distortion of the spectra 
in either | cos 0*1 or W . For overall normal- 
ization errors, fits are made with two sets of 
values of differential cross sections obtained 
by multiplying by (1 ± cre(W, | cos 0*|)), 
where a^ is the relative efficiency error; they 
are denoted as "normalizationi" in the ta- 
ble. For distortion studies, ±4% errors for 
I cos 6** I < 1 and ±4%/GeV for the W de- 
pendence are assigned, based on the un- 
certainty discussed for (9) in Sec. IV. E. 
Differential cross sections are modified by 
multiplying by (1 ± 0.08|cos6'*| =F 0.04) 
and (1 ± 0.08(W[GeV] - 1.38)) (denoted 
as "bias: I cos 6**1 ±" and "bias:VF±'', respec- 
tively). For studies of background (BG) am- 
plitudes, either hi or Ci is set to zero for ,65 
and i?i52, while either 60 or co is floated 
for Bdo- Finally, the parameters of the 
72(1270), /2(1525), and the value of r^ are 
successively varied by their errors. 

The total systematic errors are calculated 
by adding individual errors in quadrature. 
As can be seen in Table IVIIIi we obtain 



\^6(/2(1270) 



W) = 11.5 t\l tt'7 eV, 
(9) 

which is consistent with previous measure- 
ments (201. The apparent threshold en- 
hancement in the S wave is fitted in terms 
of a scalar meson, foiX) whose mass, width 
and T^~fB{rjvi) are obtained to be 

M^Y) = 1262t^78 -?03 MeV/c2 



r^„(^)=484t?^«t^«MeV 



r^^s(/o(y) 



c^ (10) 
(11) 

w) = 121 i^^r ;i^^ ev, 

(12) 



respectively. 

The mass peak of the foiX) does not coin- 
cide with the broad peak in the 5^ spectrum 



15 



in Fig. [12] due to the effects of interference. 



2. Higher mass region, up to 2.0 GeV 

Now we investigate the higher mass re- 
gion. We fix most of parameters determined 
at lower energy, and introduce, just for the 
purpose of parameterization, a single ten- 
sor resonance, f2{X), whose mass, width 
and T-yyB{rir]) are left free and fit the region 
1.16 GeV <W <2.0 GeV. We parameter- 
ize partial waves as follows: 



S 
Do 
D2 



Bdo, 



-Af„,x,e'^^^ +Bd2, 



(13) 



where Af^i^y), ^/2(i270) and ^^^(1525), are 
fixed at the values that are fitted in the low 
mass region. Here too, Bdq is set to zero 
and Bs is fixed at the values found above. 
The phases 0y , 02 and ^5 are also fixed and 
4>s ~ 0. Only the 62 and C2 parameters of 
B]j2 are floated along with the parameters 
of /2(X), i.e., its mass, width, T-yyB{rir]) and 

Two hundred sets of randomly generated 
initial parameters are prepared and fits are 
performed for each study. A unique solution 
is obtained with a fit quality of x^/ndf = 
311.4/204. The parameters obtained are 
summarized in Table IIXI Figures [13] to [15] 
show results of the nominal fit to the differ- 
ential cross sections, the total cross section, 
and spectra of 5^, Dq and Z)|. A more so- 
phisticated parameterization results in mul- 
tiple solutions. As an example, two solu- 
tions are found when the parameters of ,65 
are also floated; these are also listed in Ta- 
ble IIXI Hence we employ the simple param- 
eterization given in Eq. (13). This param- 
eterization results in discrepancies from the 
fits in some W regions for differential and 
integrated cross sections. 

Various sources of the systematic errors 
are studied and evaluated using various fits 



similar to those applied in the analysis for 
the low mass region, as summarized in Ta- 
ble [X] We take into account the errors 
for the fo{Y) parameters, as well as those 
for /2(1270) and /^(1525). We try two 
fitting regions shifted lower by two bins 
(1.08 GeV < W < 1.92 GeV) and higher 
by two bins (1.24 GeV <W< 2.08 GeV). 
For studies of background (BG) amplitudes, 
either C2 or 62 is set to zero for i?/)2 or al- 
lowed to fioat for Bdo ■ Values of cs and 65 
are changed by their errors. 

The total systematic errors are calculated 
by adding the individual errors in quadra- 
ture. The mass, width and T-y-yB{r]r]) ob- 
tained for the f2{X) meson are 



M/,(x)=1737 ±9 



+ 198 
-65 



MeV/c^ 



■h(x) 



228 



+21 +234 



T^^B{f2{X) 



(14) 
-20 -153 MeV, (15) 

+ ,,) = 5.2 tr, tTI eV, 

(16) 



respectively. 

The rather poor x^ of the fit and the 
clear disagreement in Figs. \\M and [TS] above 
the 72(1525) may imply that more than one 
tensor resonance exists in this mass region. 
Unfortunately, we cannot draw any definite 
conclusions about such a possibility from 
additional fits to the data, because inter- 
ference between amplitudes introduces too 
much additional freedom. 



VI. ANALYSIS OF THE HIGH 
ENERGY REGION ABOVE 2.4 GEV 

In this section, we present a study of the 
angular dependence of the differential cross 
section, the W dependence of the total cross 
section, the ratio of cross sections for 7777 to 
tt'^tt'^ and Xcj charmonium production in the 
high energy region, W > 2.4 GeV. 



A. Angular dependence 

As in the analysis of the tt'^tt" [11| and 
rjTi^ |12| processes, we compare the angular 
dependence of the differential cross sections 
with the function 1/ srn^O* for the data in 
the W range 2.4 GeV < W^ < 3.3 GeV. 
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TABLE IV: Parameters of the /2 (1270) and f^ (1525) assumed or fitted in Ref. [1|. 



Parameter 


h (1270) 


.f2 (1525) 


Unit 


Reference 


Mass 


1275.1 ±1.2 


1525 ±5 


MeV/c^ 


[20] 


Width 


185.lt^:« 


73t« 


MeV 


[20] 


B(f2 -> tttt) 


(84.81?:^)% 


(0.82 ±0.15)% 




[20] 


B{f2 ^ KK) 


(4.6 ± 0.4)% 


(88.7 ±2.2)% 




[20] 


B{h ^ m) 


(4.0 ±0.8) X 10"^ 


(10.4 ±2.2)% 




[20| 


e(/2 -^ 77) 


(1.64 ±0.19) X 10"^ 


(1.11 ±0.14) X 10" 


:> 


[20] 


vr 


3.62 ±0.03 


3.62 ±0.03 


(GeV/c)-i 


[2] 



5 
4.5 

4 
3.5 

3 
2.5 

t.5 

1 
0.5 



1.14 GeV 
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FIG. 10: Differential cross sections (points with error bars) and fitted curves for the nominal fit in 
the low mass region (solid curve). Dotted (dot-dashed) curves are |5p (I-D2P) from the fit. The 
vertical error bars are statistical only. 
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FIG. 11: Total cross section (points with er- 
ror bars) (|cosS*| < 1.0) and fitted curves for 
tlie nominal fit in tlie low mass region (solid 
curve). Dotted (dot-dashed) curves are \S\^ 
(|_D2p) from the fit. The vertical errors are sta- 
tistical only. 
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FIG. 12: S'^ Oo and _Di and fitted curves for 
the nominal fit in the low mass region (solid 
curve). The points with error bars are the re- 
sults of the W^-independent fits (same as those 
in Fig. 8). The vertical error bars are statistical 
only. 



TABLE V: Parameters of the /o(1370) and 
/o(1500) dOl. 



Parameter 


/o(1370) 


/o(1500) Unit 


Mass 


1200 - 1500 


1505 ± 6 MeV/c^ 


Width 


150 - 250 


109 ± 7 MeV 


B(r,r,) 


seen 


(5.1 ±0.9)% 


Z?(77) 


unknown 


unknown 



In the study of 7r°7r° data, the contribu- 
tion from the charmonia is subtracted [ll| . 
However, no rehable charnioniuni subtrac- 
tion is possible for the 7777 cross section be- 
cause of the low statistics and the larger 
charmonium component (see Sec. IV. D) 
compared to the 7r°7r° case. We limit our 
discussion in Sec. VI. A- C to the region W < 
3.3 GeV only, where the contribution of 
charmonium is small. 

Figure \W\ compares the normalized dif- 
ferential cross sections with the function, 
0.322/ sin"* 6** (solid curves). The factor in 
the numerator is calculated by dividing dif- 
ferential cross sections, which are propor- 
tional to l/sin''^^* by the total integral for 
I 008 6**1 < 0.9. Agreement is poor in the W 
region considered. A l/sin®6'* dependence 
(dashed curves in the same figure) agrees 
better with the data for W > 3.0 GeV. 
The x^'s for the 1/ sin^ 6* (1/ sin^ 9*) depen- 
dences are 29.6 (14.3) for the W = 3.05 GeV 
bin, 27.8 (7.8) for the W = 3.15 GeV bin, 
and 9.8 (4.7) for the W = 3.25 GeV bin. 
The number of degrees of freedom is 8, and 
only statistical errors are used to evaluate 
the x^- 

A l/sin^6'* dependence is not a pre- 
diction of perturbative QCD (pQCD) for 
neutral-meson pair production, and thus the 
disagreement does not imply an inconsis- 
tency with the pQCD model [1J|. How- 
ever, it nright indicate that the 7777 produc- 
tion mechanism is different from that of 
tt'^tt'^ and other production processes where 
a l/sin*6'* dependence describes data well 
for W > 3.1 GeV. The handbag model also 
predicts a l/sin''^6'* dependence for neutral 
meson pair production processes at large 
Mandelstam variable t [isl . [l6| . These pre- 
dictions are critically discussed in Ref. |27| . 
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TABLE VI: Fitted parameters for the nominal fit and for a fit without the fo{Y) in the low mass 
region. 



Parameter 



Mass (MY)) 
Width (MY)) 

r-,^B{vv){MY)) 

(f)Y 



Nominal Without fo{Y) Unit 

MeV/c^ 



1262+^g 

484l?fo 
12ll^f 







(fixed) 



MeV 
eV 
deg. 



r^7S(w)(/2(1270)) 


ii.5l^:g 


11-711;^ 


eV 


(j>2 


681^ 


66 ±4 


deg. 


-^5 


1501}^ 


164 ±13 


deg. 



bs 
cs 

62 

C2 



-2.9 
2.3 

6.5 
1.8 



+3.2 
3.5 
+2.3 
1.8 



+0.3 
0.4 



+7.6 
5.7 

+3.5 
4.6 



-8.5 
3.7 ±0.1 
-3 
6.9 



+3.6 
-4.5 
+3.0 
2.3 



^nb/GeV 

\/nb 
^nb/GeV 

Vnb 



X' (ndf) 



137.1 (119) 209.7 (123) 



TABLE VII: Three solutions obtained when rT..yS(77r;)(/2(1525)) is floated in fits to the low mass 



region. 



Parameter 



Sol.A 



Sol. B 



Sol. C 



Unit 



Mass (/o(y)) 
Width (MY)) 
T^^BivvKMY)) 



1259+^° 



1259 

502 

127 

38: 



+54 

80 
+241 
191 
139 
59 
19 
32 



1264 



+55 

86 



536_j^g3 



-162 
-69 
-18 
-30 



143 
41: 



MeV/c' 

MeV 

eV 

deg. 



r^^Z3(w)(/^(1525)) 


23.11^-^ 


8.0t?;° 


5.0+^;^ 


eV 


(j)2 


4t^° 


G8t\1 


A5tll 


deg. 


05 


188111 


1551}? 


94 ±22 


deg. 



bs 
cs 

62 

C2 



9 q+2.4 
1 r+6.0 

^^•'-'-4.8 
c 9+2.8 



-2.9 



+3.2 
3.7 



■2.3111 

3.8ttt 

3.2 ±1.2 



-2.9 
2.4 



+3.2 
3.6 
+2.6 
1.9 



^nb/GeV 



-12.5tl:l 
5.8_i;2 



^nb/GeV 



X' {ndf) 



136.4 (119) 137.2 (119) 138.6 (119) 



B. W " dependence 

We fit the W~"' dependence of the total 
cross section (| cos0*| < 0.8, where we take 
the upper boundary 0.8, to match that in 
our tt^'tt'' analysis) in the energy region 2.4- 
3.3 GeV. The fit gives 

n = 7.8±0.6(stoi)±0.4(s2/s), (17) 

and the corresponding cross section is shown 
in Fig. fTTTa) together with that of the 7r°7r° 



process in the same angular range. 

The systematic error is obtained by si- 
multaneously varying the cross section by 
±lcr at 2.45 GeV and Tlo" at 3.25 GeV, and 
by =f(M^ [GeV] - 2.85)cr/0.4 for the other 
W points in between, where a, amounting 
to 6%, is the systematic error that does 
not include the uncertainty in the energy- 
independent normalization. 

The slope parameter, ri, can be compared 
with n values in other processes that we have 
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TABLE VIII: Systematic uncertainties for the fit in tlie low mass region. 





h{Y) 


/2(1270) 


Source 


Mass 


Ttot 


r^^5(w) 


r B 




(MeV/c^) 


(MeV) 


(eV) 


(eV) 


W range 


+56.5 
-8.9 


+0.0 
-17.8 


+ 26.9 
-17.0 


+0.4 
-1.1 


BiasiW' 


+ 1.1 
-1.7 


+ 1.2 
0.0 


+ 2.3 
-2.1 


+0.1 
-0.1 


Bias:| cos6*\ 


+0.1 
-0.5 


+0.5 
-0.8 


+0.6 
-0.9 


+0.1 
-0.1 


Normalization 


+27.0 
-48.8 


+220.9 
-199.1 


+152.1 
-87.1 


+3.0 
-2.4 


BG:Bs 


+0.0 
-84.9 


+0.0 
-162.9 


+0.0 
-54.3 


+0.7 
-0.0 


BG:Do 


+49.6 
-0.0 


+0.0 
-42.1 


+ 57.3 
-0.0 


+0.0 
-2.0 


BG:I?2 


+4.6 
-10.1 


+ 100.1 
-26.7 


+30.4 
-6.0 


+ 1.9 
-0.8 


/2 mass 


+0.4 
-0.2 


+4.1 
-2.5 


+ 2.0 
-1.4 


+0.3 
-0.3 


/2 width 


+ 1.1 
-0.4 


+0.4 
-0.0 


+ 1.3 

-0.5 


+0.2 
-0.2 


/2 mass 


+8.7 
-24.7 


+41.5 
-0.0 


+21.1 
-20.3 


+2.2 
-1.4 


/^ width 


+4.8 
-5.6 


+2.5 
-8.4 


+0.5 
-2.7 


+0.3 
-0.4 


/^ T^^B{r,n°) 


+ 13.0 
-13.7 


+ 11.7 
-0.0 


+ 13.0 
-6.6 


+ 1.5 
-0.4 


TR 


+0.2 
-0.2 


+0.6 
-1.5 


+0.2 
-0.4 


+0.0 
-0.0 


Total 


+81.7 
-103.0 


+246.4 
-262.8 


+ 169.4 
-106.4 


+4.5 
-3.7 
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FIG. 13: Differential cross sections (points with error bars) in the energy bins indicated and fitted 
curves for the nominal fit in the high mass region (solid curve). Dotted (dot-dashed) curves are 
|Sp (]I52p) from the fit. The vertical error bars are statistical only. 
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TABLE IX: Fitted parameters for the nominal fit and results obtained when the Bs parameters 
are floated for the high mass region. 



Parameter 


Nominal 


Free Bs 


Unit 




Sol. A 


Sol. B 




Mass (/2(X)) 


1737 ± 9 


1742 ± 10 


1738 ± 9 


MeV/c^ 


Width (/2(X)) 


228l^J 


223^22 


236l^J 


MeV 


r,,B{vri){MX)) 


t; 9+0.9 


A 7+1-0 
*-'-0.9 


6.oll;g 


eV 


<j)x 


1591:8 


imtl 


154 ±5 


deg. 


bs 


-2.9 (fixed) 


4.4 ±0.3 


-2.8 ±0.3 


^nb/GeV 


cs 


2.3 (fixed) 
-8.1 ±0.7 


0.9 ±0.2 

q 1+0-8 


1.8 ±0.2 
-9.9 ±0.9 


^^b 


62 


^nb/GeV 


C2 


9.4 ±0.5 


9.81°;^ 


10.4 ±0.6 


v^ 


X' {ndf) 


311.4 (204) 


279.3 (202) 


288.8 (202) 


~ 



totol 


\ 




ISf 






L _._.. ID/ 






: 




++ 




/^ 


\. 


/ 






1 U 1.3 \A 


1.5 1.6 


.7 1.8 1.9 2 
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FIG. 14: Total cross sections (points with er- 
ror bars) (I cos 6* I < 1.0) and fitted curves for 
the nominal fit in the high mass region (solid 
curve). Dotted (dot-dashed) curves are \S\'^ 
{\D2\'^) from the fit. The vertical error bars are 
statistical only. 



studied earlier [j, [^, [Tl|, [T^ ■ The results are 
summarized in Table IXll The present value 
for the 7777 process is close to that for the 
^0^0 pj-Qcess, although we note the mea- 
sured W regions are different. Differences 
in this parameter among different processes 
are discussed in Ref. 12711. 



C. Cross section ratio 

The ratio of cross sections between 
neutral-pseudoscalar-meson (tt*^ or 77) pairs 



in two-photon collisions can be predicted 
relatively reliably in both pQCD and 
handbag models, based on quark charges 
and flavor-SU(3) symmetry. The pQCD 
model [IJI predictions for the cross section 
ratios for tt'^tt'^, 7777° and 7777 are summarized 
in TableEH In the table, Rf = {fn/f^of, 
where /^ (/tt) is the 77 (tt^) form factor. The 
value of Rf is not well known, and we pro- 
visionally assume it to be unity. The ra- 
tio of the cross sections is proportional to 
the square of the coherent sum of the prod- 
uct of the quark charges, |Eeie2p, in which 
ei = —62 in the present neutral-meson pro- 
duction cases. We show two predictions: a 
pure flavor-SU(3) octet state and a mixture 
with Op = —18° for the 77 and 77' mesons. 
Here, we assume that the quark-antiquark 
component of the neutral meson wave func- 
tions dominates and is much larger than the 
two-gluon component, in obtaining the rela- 
tions between the cross sections. 

The W dependence of the ratio between 
the measured cross section integrated over 
I cos 6'* I < 0.8 of 77 -^ rjri to 77 — > 
ttOttO is plotted in Fig. [Htb). For the 
tt'^tt'^ process, the contributions from char- 
monium production are subtracted using a 
mo del- dependent assumption described in 
Ref. [ll|. We use the t^t; result only below 
W < 3.3 GeV, where the charmonium con- 
tribution is negligibly small. Even though 
the ratio may have a slight W dependence, 
in order to compare with QCD (as was done 
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FIG. 15: S'^ I)o and DJ and fitted curves for 
the nominal fit in tlie iiigh mass region (solid 
curve). The points with error bars are the re- 
sults of the M^-independent fits (same as those 
in Fig. 8). The vertical error bars are statistical 
only. 



for other processes) we average the ratio of 
the cross sections over the range 2.4 GeV < 
W < 3.3 GeV and obtain 



TABLE X: Systematic uncertainties for the fit 
in the high mass region 





/2(^) 


Source 


Mass 


Ftot 


r-,-,B{rjr!) 




(MeV/c^) 


(MeV) 


(eV) 


VK-range 


+95.0 
-0.0 


+38.8 
-101.8 


+3.4 
-1.5 


Bias:W 


+3.8 
-3.8 


+14.2 
-13.0 


+0.6 
-0.5 


Bias:j cosS*| 


+0.7 
-0.8 


+6.1 
-6.1 


+0.4 
-0.3 


Normalization 


+ 125.5 
-10.4 


+ 142.7 
-0.0 


+ 16.7 
-0.0 


BG:cs 


+0.0 
-0.0 


+0.0 
-0.0 


+0.0 
-0.0 


BG:&s 


+0.0 
-0.0 


+0.0 
-0.0 


+0.0 
-0.0 


BG:Bdo 


+0.0 
-1.4 


+0.0 
-6.2 


+0.0 
-0.5 


/2 mass 


+ 1.0 

-1.1 


+4.7 
-4.8 


+0.2 
-0.2 


/2 width 


+1.1 

-1.4 


+5.1 
-4.1 


+0.2 
-0.2 


/2 r.,-,Bivv) 


+0.0 
-8.2 


+60.6 
-0.0 


+4.2 
-0.0 


4>2 


+5.1 
-9.4 


+ 11.8 
-7.5 


+0.4 
-0.0 


/2 mass 


+3.8 
-32.4 


+ 10.4 
-0.0 


+0.6 
-0.5 


/^ width 


+6.0 
-5.9 


+7.2 
-8.1 


+0.4 
-0.4 


/^ T^^Binn) 


+30.4 
-32.4 


+4.5 
-0.0 


+0.6 
-0.0 


</>5 


+9.4 
-13.3 


+24.2 
-15.2 


+0.2 
-0.0 


fo{Y) mass 


+ 14.0 
-21.8 


+93.3 
-41.1 


+6.3 
-1.9 


/(,(y) width 


+36.4 
-18.2 


+70.3 
-79.7 


+3.9 

-2.7 


U{Y) r,,B{m) 


+27.0 
-20.0 


+76.2 
-66.3 


+4.6 
-2.5 


fa{Y) phase 


+ 105.7 
-21.5 


+91.8 
-0.0 


+31.7 
-0.0 


TR 


+0.5 
-0.6 


+0.9 
-0.8 


+0.0 
-0.0 


Total 


+ 198.3 
-65.3 


+233.7 
-153.1 


+37.3 
-4.5 



']^T^\, = 0.37±0.02(stat) ± Om(sys) 
cr(7r'^7r'^) 

(18) 
for I cos 0*1 < 0.8. The prediction of this 
model with Op = -18° and i?/ = 1 
agrees well with our previous r]n'^ measure- 
ment [121, but it is in poor agreement for 
the rjrj process. However, we note that the 
W regions are different in the two cases. 

The prediction of the rjrj cross section for 
I cos 6** I < 0.6 from the handbag model is 
presented in Fig. 5 in Ref. |la|, which is 
based on measurements of other meson-pair 
production processes. We show the results 
from this measurement, which can be di- 
rectly compared with the prediction in Ta- 
ble IXIIII Agreement between the measure- 



ment and prediction is fairly good.^ 



D. Extraction of XcJ charmonium 
contribution 



As in our previous tt^^tt^ analysis |ll| . 
we extract the contributions from the XcO 
and Xc2 charmonia from the rjr] data, us- 
ing the raw yield distribution in the region 
2.8 GeV < 11^ < 3.8 GeV integrated over 



^ We do not give a quantitative comparison because 
Ref. llSlI provides only a figure without any nu- 
merical values. 
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TABLE XI; The value of n in crtot oc VK " in various reactions fitted in the W and 
indicated. The first and second errors are statistical and systematic, respectively. 



cosf^ I ranges 



Process 


n 




W range (GeV) 


1 cos 6* 1 range 


Reference 


m 


7.8 ±0.6 ±0.4 




2.4 - 3.3 


<0.8 


This work 


n.' 


10.5 ±1.2 ±0.5 




3.1 -4.1 


< 0.8 


[12] 


n\° 


8.0 ±0.5 ±0.4 


3.1 - 


- 4.1 (3.3 - 3.6 excluded) 


<0.8 


mi 


KlKl 


10.5 ±0.6 ±0.5 


2.4- 


- 4.0 (3.3 - 3.6 excluded) 


<0.6 


m 


TT + TT- 


7.9 ±0.4 ±1.5 




3.0-4.1 


< 0.6 


[3] 


K+K- 


7.3 ±0.3 ±1.5 




3.0 - 4.1 


<0.6 


[3] 



TABLE XII: Predictions and data for the cross section ratios (l4l | for tv^'tt^', rjir" and 7777 production 
processes in two-photon collisions. Here, Rf = (/ij/J^o)^, where /,, (/^) is the 77 (tt") form factor; 
the value may be taken to be Rf — 1. The rj meson is treated as a pure SU(3) octet state for the 
entries in the "octet" row, while "9p = —18°" is the most probable mixing angle between the octet 
and singlet states from experiment [20|. The first and second errors for the data are statistical and 
systematic, respectively. 



V in SU(3) 


a{n7v")/a{7T'>7v'') 


a(77r,)/a(^%«) 


Octet 
9p = -18° 


0.24i?/ 
0A6Rf 


0.367?^ 
0.627?^ 


Data (ref.) 
{W range) 


0.48 ± 0.05 ± 0.04 [12| 
(3.1 GeV < ly < 4.0 GeV) 


0.37 ± 0.02 ± 0.03 (this work) 
(2.4 GeV <W <3.3 GeV) 



TABLE XIII: Cross section integrated over 
|cosS*| < 0.6 multiplied by s^. The first and 
second errors are statistical and systematic, re- 
spectively. 



s (GeV') 


s^a(\ cos e*\ <0.6) (nbGeV^) 


6.00 


38.7 ±3.7 ±4.3 


6.50 


33.5 ±4.1 ±3.6 


7.02 


28.4 ±4.5 ±3.1 


7.56 


38.1 ±5.6 ±4.3 


8.12 


17.1 ±4.3 ±2.1 


8.70 


21.7 ±5.9 ±2.6 


9.30 


18.5 ±6.4 ±2.4 


9.92 


11.7 ±6.8 ±2.0 


10.56 


21.2 ±10.6 ±3.5 



I COS 6'* I < 0.4 (Fig. [HI), where the contribu- 
tion is enhanced against the forward peak 
from the QCD effect. 

The same formula as in our analysis for 
the 7r°7r° final state [11| is used, where par- 



tial interference between the XcO charnio- 
nium and the continuum component is taken 
into account: 



Y{W) 



{W)\^ +a{l-k)W-^, 
(19) 



+^XcJBW^,, 



where BW^^j(M^) is a Breit-Wigner 
function for the charmonium am- 
plitude, which is proportional to 



1/(W^2 



Af2 

Xc./ 



^x.j^x.j) and is 



normalized as J \BW^^j{W)\'^dW = 1. 
The masses and widths, M and F, of the 
charmonium states are fixed to the PDG 
world averages [20|. The component aW~^ 
corresponds to the contribution from the 
continuum, with a fraction k that interferes 
with the XcO amplitude with a relative 
phase angle, (p. 

We do fits with and without interference 
between the XcO and the continuum. The in- 
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terference with the Xc2 is neglected because 
of its narrow width. We assume a W res- 
olution to be O.OIW from the MC simula- 
tion, and take it into account in the fit by 
smearing the function Y{W). We apply a 
binned maximum likelihood fit with a bin 
width AM^ = 20 MeV. 

The result with interference gives nearly 
the same result as the fit without interfer- 
ence but with larger errors. The fit with 
interference cannot determine the interfer- 
ence parameters, k and </>, with a useful 
accuracy. Therefore, we take the nomi- 
nal result from the fit without interference. 
The best fit is shown in Fig. [TSl The re- 
sults are tabulated in Table IXIVl Signifi- 
cances for the charmonium signals are 5.2(t 
for the Xco and S.Ou for the Xc2- The signif- 
icances are obtained from the difference of 
the logarithmic-likelihoods with and with- 
out the corresponding charmonium contri- 
bution, where the change in the number 
of degrees-of- freedom is taken into account. 
Here, in order to obtain the most conser- 
vative value, we extracted the value in the 
interference (non-interference) case for the 
XcO (Xc2)- The systematic errors are from 
uncertainties in the W scale and the W res- 
olution (we vary them by ±3 MeV and by 
±20%, respectively) and the efficiency error. 

The results for r^-y(xcj)^(XcJ — ^ tl'n) ^-re 
consistent with the product of the known to- 
tal widths [201 and the branching fractions 
from the recent CLEO and BES measure- 
ments m, m, (8.0 ± 0.9) eV and (0.30 ± 
0.04) eV for XcO and Xc2, respectively, where 
we take the average of the CLEO and BES 
measurements. 



VII. SUMMARY AND 
CONCLUSION 

We have measured the cross section of 
77 ^ W using a high-statistics data sam- 
ple from e^e~ collisions corresponding to an 
integrated luminosity of 393 fb~^ with the 
Belle detector at the KEKB accelerator. We 
obtain results for the differential cross sec- 
tions in the center-of-mass energy {W) and 
polar angle ( |cos0*|) ranges of f.096 GeV 
(the mass threshold) < W^ < 3.8 GcV and 
up to I cos 6** I = 0.9 or 1.0, depending on W . 

The differential cross sections are fitted 




0.2 0.4 OS 0. 



|cos 9*1 

FIG. 16: The angular dependence of the differ- 
ential cross sections in different W regions, with 
the normalization to the cross section integrated 
over I cos9*\ < 0.9. The solid and dashed curves 
are proportional to l/sin^6* and l/sin^O*, re- 
spectively, normalized similarly. 



in the energy regions, 1.12 GeV< W < 
1.64 GeV and 1.20 GeV< W < 2.00 GeV 
using a simple parameterization of the S, 
Dp and D2 waves, assuming that ampli- 
tudes consist of resonances and a smooth 
background. In the low energy fit, con- 
sistency of the parameters of the /2(1270) 
and /2(1525) with previous measurements 
is checked. The apparent threshold en- 
hancement in the S wave is fitted in terms 
of a scalar meson, fo{Y), whose mass, 
width and T^~/B{rir]) are obtained to be 
1262 +Jrl tll^ MeV/c^ 484 tl^o tjtl MeV 
and 121 tsf till eV, respectively. The 
foiY) is introduced only to parameterize the 
data and may not be a single resonance. 

For the energy region of 1.20 GeV < 
W < 2.00 GcV, fits are then performed 
by fixing most of the parameters obtained 
in the low energy region and by including 
an additional tensor resonance. The ob- 
tained mass, width and T^^Bijiri) for the 
tensor meson are 1737 ± 9 tgf MeV/c^, 
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TABLE XIV: Charmoniuin yields and r^.yB{riri) from the present measurement. Two cases are 
shown: with and without interference between XcO and continuum. The first and second (if given) 
errors are statistical and systematic, respectively. Only differences in log-likelihood values are 
meaningful. 



Interference 



Yield(xco) Yield(xc2; 



-2 In L/ndf 



ri-t{Xco)t3{xco 



mi) (eV) r^^(xc2)Z3(xc2 -> w) 



Without 
With 



21.7 ±5.3 8.5 ±3.6 
21.5 ±9.2 10.1 ±3.9 



39.5/46 
38.5/44 



1.4 ±2.3 ±1.2 



0.53 ±0.22 ±0.09 




0.8 
pO.6 
lo.4 

0.2 




(b) 



-^-r 



2.4 2.5 2.6 2.7 2.8 2.9 3 
W(GeV) 



FIG. 17: (a) The W dependence of the cross 
sections (\ cos9*\ < 0.8) for the tt^tt" (open 
squares) (lll | and Tjri (closed circles) processes. 
The curve is the power-law fit for the latter pro- 
cess, (b) The W dependence of the cross section 
ratio of r]ri to 7r°7r° (| cos 6**1 < 0.8). The line is 
the average in the 2.4 - 3.3 GeV range. The er- 
ror bars are only statistical in the above figures. 




3 3.2 3.4 

W(GeV) 

FIG. 18: The W distribution in the charmo- 
nium region (|cos0*| < 0.4). The fit does not 
take interference into account. 



228 tfo -lit MeV and 5.2 +°;^ tlY eV, re- 
spectively. The f2{X) is a parameterization 
used to describe the data in 1700 MeV mass 
region. It may represent some of the possi- 
ble tensor resonances in this mass region. 

We observe clear signals from /2(1270) -^ 
■qrj and /2(1525) — > 7777 for the first 
time in two-photon collisions. The 

product T^^B{rjrj) for the /2(1270) is 
11.5 Is'o ilr eV. Our J2{X) may corre- 
spond to the /2(1810) state reported in 
Ref. [2J|. The result of our measurements 
for the product T-y.,B{'qri) for the /2(1270) 
and /2(1525) are consistent with the previ- 
ously known values [20 - 24| - 

The angular dependence of the differen- 
tial cross section in the 2.4-3.3 GeV region 
are compared with ~ 1/ sin 6* dependence, 
as found in the tt^tt" process ll l| and pre- 
dicted by the handbag model jlSl . Il6| for 
VT > 3.1 GeV. However, in the 7777 process, 
a 1/ sin^ 9* dependence is not found in the 
data for the the energy region where the 
measurement is performed. 

The slope parameter n for the cross sec- 
tion, a{W) ^ W~"', in a similar W region 
is close to that measured in the tt^tt" pro- 
cess (ll| . 

The measured cross section ratio, 
a{r]r])/a{TT°TT°) = 0.37 ± 0.02 ± 0.03 (for 
I cos 6**1 < 0.8), is compared with the pre- 
diction of pQCD [IJI with a pseudoscalar 
meson mixing angle, 6p = —18°. We 
find that the assumption for the squared 
form-factor ratio, Rf = {fri/f-K°Y — 1j 
which is in a good agreement with the ratio 
air^-K^^lairP-K^) [iJl cannot reproduce well 
the ?7?7 measurement. Our result agrees 
rather well with the recent handbag model 
prediction |16| . 

Gharmonium contributions in the 7777 
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process are confirmed for the first time. 
Our measurements are consistent with the 
known partial decay widths of the XcO and 
Xc2 to 77 [201 and 7777 [2^, [2^ final states. 
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